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Summary 

A code for computing the aerodynamic design of a 
multistage axial-How compressor and, if desired, the 
associated blading geometry input for internal flow 
analysis codes is presented. The aerodynamic 
solution gives velocity diagrams on selected 
streamlines of revolution at the blade row edges. 
Blading is defined from stacked blade elements 
associated with the selected streamlines. The blade 
clement inlet and outlet angles are established 
through empirical incidence and deviation angle 
adjustments to the relative flow angles of the velocity 
diagrams. The blade element centerline is composed 
of two segments tangentially joined at a transition 
point. The local blade angle variation of each 
segment can be specified with a fourth-degree 
polynomial function of path distance. Blade element 
thickness also can be specified with fourth-degree 
polynomial functions of path distance from the 
maximum thickness point. 

Steady axisymmetric flow js assumed; so the 
aerodynamic problem can be reduced to solving the 
two-dimensional flow field in the meridional plane. 
Because the equations of motion as developed herein 
are only applicable for calculation st'Mions outside 
the blade rows, stations at the blade edges, but not 
inside the blade rows, are used. The streamline 
curvature method is used for the iterative 
aerodynamic solution. If a blade design is desired, 
the blade elements are defined and stacked within the 
aerodynamic solution iteration. Thus the design 
velocity diagrams can be located at the blade edges. 

The program input includes the annulus profile, 
the overall compressor mass flow, the pressure ratio, 
and the rotative speed. A number of parameters are 
input to specify and control the blade row 
aerodynamics and geometry. There arc numerous 
options controlling the way information Is input 
and for specifying the amount of output. The output 
from the aerodynamic solution has an overall blade 
row and compressor performance summary followed 
by blade element parameters for the individual blade 
rows. If desired, blade coordinates in the streamwise 
direction for internal flow analysis codes and/or 
coordinates on plane sections through blades tor 
fabrication drawings can be printed and punched. 


Introduction 

The axial flow compressor is used for aircraft 
engines because it has distinct configuration and 


performance advantages over other compressor 
types, bin the good potential performance is not 
easily attained The problem and challenge to the 
designer is to model the actual flows well enough to 
adequately predict aerodynamic performance. 
Progress is continually being made with codes for 
computing the complex three-dimensional flows in 
turbomachinery. However, it is extremely difficult to 
design mechanically acceptable lurbomachinery 
blading by using the direct approach (i.e., specifying 
inviscid blade surface velocities and computing the 
blade geometry). Consequently, the more detailed 
codes are generally used in the analysis mode; that is, 
the flow field is calculated for a fixed geometry. The 
current procedure is to establish blading geometry 
with simpler design codes and then to use the more 
detailed analysis codes in blade rows where 
troublesome flow conditions are likely to exist. In 
this way prototype designs can often be adjusted 
before hardware is built and tested. 

The time and effort needed to get acceptable 
configurations can be reduced if the design code can 
be made to yield a good initial solution and if the 
design and analysis codes can be made more 
compatible with one another. This compatibility can 
be achieved (1) if the output from a design code can 
be directly used by flow and mechanical analysis 
codes and (2) if corrective adjustments indicated by 
the analysis codes can effectively be made in the 
design code. With these objectives in mind a 
composite aerodynamic and blade design code for 
axial-flow compressors has been developed. The code 
and its capabilities are the subjects of this report. 

The aerodynamic solution assumes steady, 
axisymmetric flow and uses a streamline curvature 
method for calculation stations outside the blade 
rows. The program is structured so that the empirical 
correlations (such as those for loss, deviation angle, 
and incidence angle) can readily be changed when the 
need or desire exists. The method of descrioing 
blading is a compromise between the vast amount of 
input needed for completely general blade elements 
and the restrictions of simple shapes. .A blade 
element is defined on a conic surface with thickness 
applied to a centerline that is composed of two 
segments tangentially joined at a transition point. 
The blade angle function of each segment can be 
defined with a fourth-degree polynomial. Thickness 
is prescribed by first specifying a maximum thickness 
value and location. The distribution of thickness in 
each direction from the maximum thickness location 
is then prosciibcd with a tourth degree polynomial, 
f inally each polynomial coefficient is defined across 







blade dements with a third-degree polynomial 
function of annulus height. 

Compressor Design Procedures 

The discussion of the compresso*- design 
procedures is organized according to usage in the 
computer program; so for better orientation an 
operational overview of the program is given first 
(table I). The computer program can be divided into 
three major phases of calculation: (I) the input and 
initialization phase. (2) the iteration phase, and (3) 
the terminal calculation phase. In the input and 
initialization phase the input data are read and 
interpreted, the calculation stations are located with 
estimated values for the blade edges, and streamlines 
are located on the basis of annulus area. Estimates of 
stagnation temperature and pressure and axial and 
tangential velocity components are also made for all 
calculation points in the flow field. 

The iteration phase includes both the flow field 
and the blade design iterations. In the flow field 
iteration the equations of motion are satisfied in the 
meridional (r-c) plane for stations that are lines 
across the flow annulus. At the stations the equations 
of motion and overall flow continuity arc satisfied 
with fixed values of streamline slope and curvature 
for a complete computational pass across the 
annulus. After the overall flow continuity condition 
at a calculation station is satisfied, the internal 
streamline intersections with the station lines are 
updated by solving for the locations that give 
specified fractions of overall station weight How. At 
the completion of a pass through all the calculation 
stations in the annulus, the new streamline locations 
are curve fit for new streamline slope and curvature 
values. 

To insure proper location of the blade edge 
stations, most of the blade design iteration is made 
concurrent with the flow field iteration. This 
operation includes the calculation of incidence and 


deviation angles, the layout and stacking of blade 
elements, and the realignment of the elements. 

The terminal calculation phase performs the final 
calculations and generates the output . Mass averaged 
parameters for the individual and cumulative 
compressor blade rows are computed and printed 
first. Then tabulated values of aerodynamic and 
blading parameters along the station lines are 
computed and printed Finally blade section 
cix)rdinates and other sectioii mechanical properties 
can be computed and printed if desired. 

The program is discussed in greater detail in the 
following subsections. 

Input Initializnlion 

The basic computational plane is the meriodional 
(r-z) plane of a cylindrical coordinate system. A 
graphic view of an example compressor 
configuration is shown in figure 1. The hub and tip 
casing walls are fixed input. Calculation stations arc 
located at the blade row leading and trailing edges 
and at other annular locations for the purpose of 
locating streamlines. The input data can be classified 
into two groups: general information and calculation 
station and blade row information. The input 
parameters and options, along with the input data 
format, are described in appendix B. (All 
mathematical symbols are defined in appendix .A.) 
Additional advice on how to set up the input is given 
in the section User Information. 

Oneral Information 

All the general information is read in first. 
Included are the following: 

(1) Compressor rotational speed 

(2) Inlet flow rate 

(3) Desired compressor pressure ratio 

(4) Gas molecular weight 

(5) Number of streamlines 
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(6) Number of biade rows 

(7) Number of annular stations 

(8) Coefficients for Cp as c fifth-degree 
polynomial function of temperature 

(9) Far upstream values of total temperature, 
totad pressure, and inlet tangential velocity for each 
streamline 

(10) Streamtubc mass flow fractions between 
streamlines 

(11) Sets of points to define tip and hub casing 
contours 

(12) Sets of blade element profile loss parameters 
that are tabulated as functions of blade element 
loading parameter and fraction of passage height 

As many as five loss sets can be input. The particular 
loss set used for a given blade row is designated in the 
blade row input. Usually at least two loss sets are 
input — one for rotors and another for stators. 

CidcidatkHi Statioa Data Sets 

The data sets that contain information about the 
calculation stations and blade rows arc read in order 
from annulus inlet to outlet. The first card of the 
data set identifies the type of station, the tip and hub 
axial locations, the tip and hub boundary layer 
blockage factors and the station mass flow bleed. For 
annular stations the single card is the whole data s>n . 
For rotors and stators several cards are used to 
describe (1) the blade row inlet and outlet station 
information, (2) the blade row aerodynamic 
parameter input and controls, and (3) the parameters 
defining blade geometry. A blade and the associated 
edge calculation stations are located in the annulus 
by using a reference blade element stacking line. 
Stacking axis tip and hub axial locations and lean 
angle in the circumferential direction are input. 

The locations of the calculation stations at the 
blade edges are at first approximated from some of 
the input blade geometry information. The station 
locations are moved during later iterations when the 
blade elements are defined and stacked. However, 
the input tip and hub boundary layer blockages and 
mass flow bleeds for the inlet and outlet stations are 
constant. 

Aerody iamic paramettr input and controls . — The 
blade aenxlynamic design is controlled with several 
parameters that impose the necessary and sufficient 
conditions for a solution. The options as to how such 
conditions can be imposed are shown in table 11. For 
rotors the most convenient option is to specify the 
stage energy addition as a cumulative fraction of the 
overall compressor energy addition. With this option 
the radial distribution of energy addition is not input 
directly but is imposed throu^ a normalized rotor 
exit stagnation pressure profile that is expressed as a 


polynomial function of annulus height in the radial 
direction. The pressure level is computed internally 
to the program from the energy input level and the 
computed losses. With the other rotor options the 
exit temperature profile is input instead of the energy 
addition fraction being specified. For either a rotor 
or stator, stagnation pressure profiles can be input 
instead of the losses being computed internal to the 
program. These options can be useful to users who 
have existing aerodynamic designs but want to use 
this program for blade description and fabrication 
coordinates. 

At a stator exit a tangential velocity profile is input 
as a polynomial function of radius. Unless specified, 
the stator outlet pressure profile is determined from 
stator losses and streamline mixing effects from the 
upstream station. 

There are some input aerodynamic limits that the 
program will not allow t be exw ' ^. For a rotor 
the limiting parameter; e tip difl-.».on factor and 
absolute flow angle jx the hub. The stator 
aerodynamic limits are diffusion factor and into 
Mach number at the hub. If an aerodynamic limit is 
exceeded during iteration, the stage energy addition 
is lowered by the amount needed to get the 
aerodynamic limit within bounds. If any other stage 
is not up to one of the aerodynamic limits, the energy 
decrement is made up among such stages. If all the 
stages reach an aerodynaruc limit, the input overall 
compressor pressure ratio is lowered. 

The blade angles are related to fluid flow an^es 
along streamlines by two key correction 
parameters — incidence angle at the into and 
deviation angle at the outlet (fig. 2). 1 here are several 
options for specifying each. Two of the options for 
both the incidence and deviation angles are the two- 
and three-dimensional methods of reference 1. Tly 
other incidence angle option is user-entered tabular 
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data referenced to either the centerline or the suction- 
surface blade angles at the inlet. Other deviation 
options are user-entered tabular data and a version of 
Carter’s rule, which was modified to account for 
centerline shapes other than a circular arc. The 
modification is shown in figure 3. 

Another input aerodynamic parameter is the 
minimum blade choke margin (/I//4*) - I , where A is 
the local streamtube cascade channel area and ^4* is 
the corresponding area needed for choked flow. The 
A* value is the area needed to pass the streamtube 
flow at a relative Mach number of i.O. The effects of 
losses in all blade rows and energy addition in rotors 
are included in the computation of A*. Choke 
margin depends on the flow conditions and geometry 
defining the channel area. If insufficient choke 
margin exists in a prototype design, some 
compromise must be made in either the aerodynamic 
requirements or the geometry. Minor choke margin 
deficiencies can usually be accommodated with 
adjustments in geometry. Logical procedures for 
geometry adjustments are not obvious; however, if 
the minimum margin occurs at the channel entrance, 
increased incidence is an effective method of relief. If 
a minimum desired choke margin is input, the 
program will adjust incidence angle up to +2“ to the 
leading-edge suction surface in order to attain the 
specified choke margin if the channel entrance is the 
problem. When the minimum margin occurs at other 
locations in the channel, the minimum value and its 
location are printed in the output and it is up to the 
user to decide if he wants to make compromises to 
improve the choke margin. 


Blade geometry parameters. — A number of blade 
geometry parameters are input for the pupose of 
defining a blade. Blade chord is defined ^ong flow 
streamlines, but for the purpose of this blade 
definition a radial projection of streamline chord is 
specified because it is more meaningful for defining a 
structurally sound configuration. The radially 
projected chord is defined from the number of 
blades, the tip solidity, and a normalized polynomial 
for the radial variation of chord. The blade is 
basically defined from a stacked series of gradually 
changing airfoil shapes or “blade elements” in the 
radial direction. 

Each blade clement, as shown in figure 4, is 
defined from a thickness distribution applied to a 
two-segment centerline. The variation of the local 
centerline angles #c with path distance can be specified 
by option through the parameter IDEF(IROW). If 
IDEF(IROW) equals zero, the k for each segment 
varies linearly with path distance (as a circular arc). 
When IDEF(IROW) docs not equal zero, the k for 
each segment is expressed as a fourth-degree 
polynomial function of path distance. The blade 
angle is continuous at the transition point, but the 
rate at which the angle changes with distance 
(curvature) can be discontinuous. The ratio of 
curvature for the first segment to that for the second 
segment is defined as the turning rate ratio. When the 
blade local centerline angle k is specified by 
polynomial coefficients, the turning rate ratio is 
controlled by the relative magnitudes of the linear 
term coefficients of the polynomials for each 
segment. However, when the segments arc treated 
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simply as circular arcs, the turning rate ratio is a 
blade element input parameter. 

When IDEF(IROW) equals zero, there are some 
options for specifying the turning rate ratio at the 
transition point. With the CIRCULAR option the 
value is set at 1 .0. as for a circular arc blade element, 
for all blade elements in the blade row. With the 
TABULAR option a table of values for the elements 
is read. With the OPTIMUM option a value will be 
set by an empirical function of inlet relative Mach 
number. For this option the blade element will be a 
circular arc below a relative Mach number of 0.8. As 
relative Mach number increases, the ratio of first- to 
second-segment turning rate at the transition point is 
reduced. A limit of zero camber on the suction 
surface of the first segment is app. oachcd at an inlet 
relative Mach number of about 1.60. 

The coefficients for the centerline polynomial (i.e., 
when IDEF(IROW) 0) are input as a cubic funrtion 
of blade span. There are two reasons for this method 
of specification. First, the user is more confident of 
specifying a relatively smooth blade surface; and 
second, the amount of input is reduced over that 
required by individual coefficients for as many as 1 1 
blade elements. 

Blade element surface definition begins with three 
anchor points from the centerline. These points are a 
maximum thickness point and the two end points. A 
maximum thickness value normalized by chord and 
its location as a fraction of chord are input. At the 
maximum thickness point the normal-to-centerline 
distance to each surface is one-half the maximum 
thickness, and the surface « angles are equal to the 

centerline k. . 

At the blade element ends the leading- and trailing- 
edge end circle radii normalized to chord are input. If 
IDEF(IROW) does not equal zero, the end 
configurations are ellipses with semimajor axes 
tangent to the local centerline. For this case the input 
end circle radius is used as the minimum radius value 
of the ellipse. For each ellipse c ne other parameter is 
input to speci^'y elongation. The parameter is 
e = (b/a)-\, wi'ere b and a are the semimajor and 
semiminor axes, respectively. Note that as e 
approaches zero, the ellipse approaches a circle with 
the input radius. 

A surface definition criterion is that the surface 
curve join the end circles or ellipses at a point of 
tangency. When IDEF(IROW) equals zero, the 
surface curves are defined with k being a linear 
function of path distance for each segment. As 
explained in reference 2, necessary and sufficient 
conditions exist to completely define the surfaces 
when the computation is begun on the segment where 
the maximum thickness occurs. 

When lOEF(lROW) does not equal zero, the blade 


surfaces for each segment are defined by polynomial 
distributions of the normal-to-centcrline distance. 
The functional relation for this distance is 


r= -y -a\ S„ -t-a\S- Sf,+ 


-bS~ -cS^-dS* 


where S is the centerline distance (normalizeu by 
chord) from the maximum thickness point. Values of 
S are positive in either direction from the maximum 
thickness point; and Sq is the maximum S, which is 
the distance to the point where the end ellipse 
intersects the centerline (fig. 4). 

There are two other input parameters for blade 
rows. One is a material density for rotors. If a 
nonzero value is input for a roto.*, the stacked blade 
will lean in both the meridional and the r-6 plaiies so 
that the centrifugal force on a blade with the input 
material density will balance the aerodyn^ic forces 
at the design point. The objective is to minimize the 
blade root stress. With atmospheric air as the 
working fluid, the lean is normally only a fraction of 
a degree. 

The final input parameter, NXCUT, controls the 
number and location of planes through a blade row 
for which fabrication coordinates are desired. If the 
parameter is zero, the program will set the number of 
XCUT’s on the basis of aspect ratio, which is the 
ratio of overall radial to axial blade lengths. For 
positive parameter input values the program will 
determine appropriate locations for that number of 
planes to represent the blade. Negative parameter 
values trigger an option to read cards for the XCUT 
plane values. The number of input values expected 
for a blade row is the absolute value of the negative 
parameter. 

(nitialization 

Once the input is read, a number of initialization 
calculations are made in subroutine START in 
preparation for the iterative phase of computation. 
The axial locations of the blade edges are 
approximated and the intersections of all station lines 
with the casing walls are determined. Checks ^e 
made to be certain that the spacing of calculation 
stations is appropriate. Annular stations will be 
shifted by the program if calculation stations cross 
one another or if adjacent spacing is less than 30 
percent of the spacing of neighboring stations. 

Streamlines are initially positioned by applying 
the input stream-tube weight flow fractions to the 
annulus area. From the input data the cir- 
cumferential component of velocity and the 
stagnation temperature and pressure are 
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foot at the hlade surface because the boundary layer 
cannot sustain the sudden static pressure rise. In 
either case the shock losses are less than those 
predicted by a normal shock. To empirically account 
for these effects, the computed normal shock loss is 
reduced by dividing by the average iniet relative 
Mach number squared. 

All the other blade row losses— profile, secondary, 
etc.— are represented by a correlation with fraction 
of passage height and aerodynamic blade loading. 
The values for such a correlation are input in tabular 
form. The aerodynamic blade loading parameter in 
the table is the diffusion factor of reference 1 . In 
equation form it is 
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The loss parameter in the table is 
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where •-> is the loss coefficient. 
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The rotor e.vit tangential velocity is calculated 
directly from the Euler equation 


A velocity gradient procedure is used to construct 
the axial velocity profile from the tip to the hub with 
the stagnation state values, the streamline 
characteristics, and the tangential component of 
velocity held fixed. Since this inner loop of the 
program is used many times, some effort was made 
to evaluate its accuracy and efficiency for typical 
streamline spacing. Reasonably good accuracy and 
stability were found to result from a rather simple 
procedure. Let 

( 8 ) 
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H2 -Hi = Cp dt = U2\9, - ^ ( 5 ) 

*■ 

Note that the enthalpy change is evaluated by using 
an integral for the calorically nonperfect gas; that is. 
Cp is a function of temperature. All state processes in 
the program use thermally perfect, but calorically 
nonperfect, gas relations; so integrations and in some 
cases iterations are used in several small function 
routines. 

Inner loop.— The basis function of the inner loop 
is to determine the axial velocity profile at the 
calculation station. The axial velocity level is set by 
flow continuity, and the distribution is controlled by 
the radial equation of motion. The differential 
equation 'S developed in appendix C. The form used 
in the program is 
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With a and b constants for the / interval along the 
station path, the solution for is 
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A two-step procedure is used in the program. First 
a, b, and V„ values on the streamline j are used to 
aetermine a temporary i . The a and b values 
are slightly dependent on SO ^mj + 1 is used to 
determine new a and b values. The second step uses 
the average of the old and new respective values of a 
and b to compute a final 1 value. This Vmj^x 
aiue will then be used as the current value for 
*e next / interval. 

When values are set on ail streamlines, flow 
continuity is checked by using dr mtegration of a 
piecewise cubic curve fit of pV^r values at the 
streamlines. If the integrated weight flow is not 
within 0.01 percent of its specitied value, the tip 
reference Vm is adjusted and the V„ profile is 
reconstructed. The method of adjusting the reference 
value of is shown graphically in figure 5. There 
arc iwo solutions to the continuity equation in 
compressible flow — the subsonic and supersonic 
solutions. When a parabolic lit of trial solutions is 
used to get a new trial value of the lower or 
subsonic solution is always sought. The V„ 
adjustment between iterations usually is small; so 
convergence normally is achieved in three or four 
passes. 

Once convergence is achieved, the profile is back 
integrated to find the fraction of wei^t flow points 
represented by the streamlines. These points are 
saved until the outer loop pass through all the 
stations is completed for the purpose of relocating 
streamlines. 

Blade Design 

A blade is defined from stacked blade elements. 
The procedure for laying out blade elements and 
stacking them for blade definition is given in detail in 



reference 2. Only a summary description is given 
herein. A blade element is laid out on a cone with a 
center axis coincident with the turbomachine axis of 
rotation. The angle and location of the cone are fixed 
by the intersection of the streamline with the leading* 
and trailing-edge station lines of the blade (fig. 6). 

The leading- and trailing-edge blaue angles are 
related to aerodynamic flow angles primarily through 
two key correlation parameters — incidence angle and 
de'.’idtion .ingle. The user has some options for the 
specification of these correlation parameters, as 
already discussed in the section on data input. 
Application of incidence and deviation angles to the 
flow angles at the blade edges gives blade angles in 
the local streamwise direction. Corrections to 
“cascade” deviation angle for a change in radius and 
axial velocity are made internally to the program. 
These corrections are presented in reference 4 to 
relate deviation angle to a cascade section with 
equivalent circulation rather than with the same 
camber angle. 

Because the cone angle of the associated blade 
element is usually a little different from these local 
streamwise blade angles, corrections are made with 
current streamwise and radial direaion derivatives. 
The blade element reading- and trailing-edge angles 
are calculated frem aerodynamic flow angles in 
subroutine BLADE. 

Blade element layout . — There are several options 
for controlling the blade element layout (see the 
IDEF (IROW) parameter description in appendix B). 
With all but one of these options a blade element is 
described by a prescribed thickness applied to a 
prescribed centerline (fig. 4). The centerline is treated 
as two segments that are joined at the reference 
transition point. The rate of change of the local blade 
aiigk with path distance, n-fis) (fig. 4), is controlled 
by a fourth-degree polynomial for each segment. The 
coefficients for the polynomials are input, but they 
are scaled in the program to match blade element 
inlet and outlet angles. The fourth-degree polynomial 
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representation of segment blade angle represents 
greater specification freedom than does the linear 
specification of reference 2, wliere the ratio of inlet- 
to-outlet segment curvature at the transition point is 
input rather than any polynomial coefficients. A 
summary derivation of the equations for the 
centerline coordinates is given in appendix D. 

Blade element thickness is defined along a path 
that is locally normal to the centerline. The pressure 
and suction surfaces are equidistant from the 
centerline. Thickness is specified in both the forward 
and rearward directions from the maximum 
thickness point by polynomials of the form 






( 12 ) 


The input coefficient': are scaled to meet the hading- 
and trailing-edge ellipses at the appropriate tangency 
points. The control routine for the blade element 
layout in uie program is CONIC. 

B/ade element stacking,— Iht rotating parts of 
turbomachinery normally operate at high stress levels 
because of high centrifugal force. The high 
centrifugal acceleration also causes stress from 
bending moments to be very sensitive to blade 
element location. Thus it behooves the designer, 
first, to be reasonably accurate m the stacking 
compulation and, second, to try to minimize stresses 
that can be easily reduced — namely, those from the 
steady-state bending. The blade bending moments 
from aerodynamic forces can be counterbalanced by 
centrifugal force moments with slight blade lean in 
both the (r-z) and (r-0) planes. 

The reference line for stacking purposes is a radial 
line through the hub stacking reference point (fig. 7). 
The sections used for stacking alignment are planes 
normal to this reference line in space. Such planes are 
used because their centers of area are essentially the 
centers of centrifugal force also. The slacking line is 
a line that can be leaned from the reference line at the 
hub reference point. For alignment purposes the 
planes pass through the stacking line intersection of 
blade elements (fig. 0- Blade sections are defined by 
interpolation across blade elements. When the 
section center of area does not match the stacking 
line, the corresponding blade element is translated 
and rotated on its cone for the stacking adjustment. 
Normally the adju:>tmenis decrease by about an order 
of magnitude for successive passes through the 
stacking procedure. For each pass the stacking axis 
lean angles in both the (r-z) and (r-6) planes are 



Figure 7. - location of blade sections for blade element 
stacking adjustments. 


recomputed and adjusted if the stacking axis lean 
option is activated through the input data. 

Terminal Calculations and Output 

The program output of an example two-stage 
compressor is shown in table III. In general the 
output is printed shortly after its computation so that 
large arrays of data are not stored. Data are printed 
from each of the major phases of 
computation— input, iteration, and terminal. The 
first information (tabic III (a)) is the input data, 
which are printed directly from input routines in very 
nearly the order in which the input was read. 

The second major part of output (table HI (b)), 
from the iterative phase of comt utation, is printed to 
help the user monitor the solution. Although these 
data have little value once the solution is converged, 
they are quite helpful in disclosing bad input and in 
finding sources of problems when solutions arc not 
achieved. 

For computational stability a station aspect ratio, 
defined as (rf -Z/). is limited to 7 for 
streamline fits. When the limit is exceeded, particular 
stations (according to the priorities set forth in the 
section User Information) are eliminated from the 
curve fils used to locale streamlines. The first data 
shown from the iteration phase are a table ot such 
calculation station infornnation (table II (b)). On the 


left is a list of calculation station locations used to 
compute streamlines, along with the «::ociated 
aspect ratios. On the right is the input list of station 
locations and aspect ratios. When blade rows are 
Slacked, tiie blade edge stations are relocated, and 
thus the station aspect ratios change. After the first 
stacking on iteration 4, the station aspect rat: is are 
rechecked and changes in the station list are made if 
necessary. 

Arrays of axial velocities throughout the flow field 
for each iteration are the bulk of the output printed 
from the iterative phase of computation. These data 
are useful for observing solution stability since the 
solution convergence criterion is based on changes of 
axial velocity between successive iterations. Some 
compressor overall parameters are shown above the 
velocity arrays. Parameters included are the overall 
values of input pressure (PR), current computed 
pressure ratio (CPR), enthalpy increase (DHC), and 
ideal enthalpy increase (DHI). 

When the aerodynamic solution is converged, the 
overall parameters for individual blade rows and the 
overall cumulative values in the compressor are 
computed and printed. Overall temperature and 
pressure values are calculated by mass averaging their 
equivalent enthalpy values. The cumulative forward 
axial thrust is the axial force exerted on the rotating 
shaft by aerodynamic forces from the hub inlet 
station of the first blade row to the local point. The 
thrust force shown for individual blade rows is the 
axial force on the shaft from the trailing edge of the 
upstream blade to the trailing edge of this blade row. 
3mce the blade forces on stationary blade rows act on 
the casing, the thrust value on the rotating shaft is 
simply the static pressure force on the tapered shaft 
in the forward axial direction. Effects of cavities 
below the hub flow path are not included since 
undetermined information about seal locations and 
pressure differences would be needed. The gas 
bending moments are values for a single blade. The 
bending moments are referenced to the stacking axis 
intersection with the flow path wall from which the 
blades are attached. 

Sets of calculation station data for streamlines 
across the channel follow the overall data. For all 
stations, velocity components, streamline slope and 
curvature, and both stagnation and static values of 
temperature and pressure are given. For stations at 
blade row edges, additional information is computed 
and printed. These parameters are (Da complete 
description of velocity triangles, (2) definition of 
blade eiemen»s, (3) relations between aerodynamic 
and blade angles, (4) aerodynamic performance 
parameters, (5) streamline choke area margin, (6) 
local blade force intensity in pounds per radial inch 
on a blade, and (7) blade edge direction derivatives 
r dBfdr 



Figure 8. - Coordinate system for blade section | 

output data. | 


If the input options call for fabrication 
coordinates, they are printed after all the 
aerodynamic output. The coordinates are printed in 
tabular form with four sections on a page, as shown 
in table 111(c), The length coordinate L is a distance 
along the chord line, with the most forward point 
being zero (fig. 8). The pressure- and suction -surface 
height values Hp and respectively, are referenced 
from the chord line. Surface height values are given 
for at least 20 round-value increments of L; also 
surface coordinates are given for three specific values 
of L— the blade trailing edge and the leading- and 
trailing-edge ellipse tangency points with the 
surfaces. 

A blade section’s properties are shown above its 
table of coordinates (table 111(c)). The blade section 
radial location, the L and H stacking point values, 
and the section setting angle are given to locate and 
orient the blade section. The blade section center-of- 
area coordinates, section area, minimum and 
maximum moments of inertia through the center 
area, orientation angle of the maximum moment of 
inertia with respect to the axial direction, section 
torsion constant, and twist stiffness are all useful 
information for design and stress analysis. 

After all the fabrication coordinates for a given 
blade row arc printed, the blade section coordinates 
are presented in another orientation that may be 
more useful for further flow analysis. With a 
stacking axis reference, coordinates for the same 
blade sections are given in the axial and tangential 
directions. 


User Information 

Since earlier sections of the report discuss the 
input, output, and main centers of program control, 
this discu.>sion is directed at the user who is trying to 
get the program on his computer and to make it run 
efficiently. Some facts about the program as well as 
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some advice about the input are given. 

The code, which is written in FORTR.V.N, takes 
about 80,000 decimal words of computer storage. 
The call relation among the “subroutines is shown 
graphically in figure 9. Note that the tickmarks on 
the routine boxes in the figure mean that there are 
other call lines to the routine. These lines are shown 
on the other part of the figure where the routine 
name is repeated. The program running time on 
either i Univac 1110 or an IBM 360-67 is about 2 
minutes for a single-stage compressor and about 5 
minutes for a five-stage compressor. Several of the 
key indices in COMMON/SCALAR/ are described 
in the following tabulation. 

Index Description 

1 calculation station index after 

preliminary calculations are 
completed. The program is 
dimensioned for 50 calculation stations 
and 20 blade rows, of which only 10 
can be rotors. Each blade row accounts 
for two calculation stations — one at 
the leading edge of the blade and the 
other at the trailing edge. Rotors, 
stators, and annular calculation 
stations can be put together in any 
combination with the following 
constraints: The number of stations 
cannot exceed 50. There must be at 
least four annular stations ahead of the 
first blade row and at least three 
annular stations behind the last blade 
row. 

IROTOR rotor index 

•ROW blade row index 

J streamline index. Streamlines are 

numbered from one at the tip. 

K loss set index for subroutine INPUT 

.\s indicated in the table at least four annular stations 
are expected upstream of the first blade row and at 
least three downstream of the last blade row. 
Additional annular stations can be located between 
blade rows but not within blade rows; that is, not 
between the inlet and outlet stations of a given blade 
row. 

Streamline intersections of station lines are 
determined by integrating velocity profiles at station 
lines to the specified ma.ss flow fraitions. Streamline 
slope and curvature are determined from streamwise 
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curve fits of these intersections. The consequence of 
this procedure is that the number of iterations .nd 
the program convergence characteristics are 
dependent on the calculation station location 
although the final solution, in general, is not very 
dependent on the location of the calculation stations. 

The user can reduce the number of iterations and 
hence the program running time with good placement 
of calculation stations. The first calculation station 
should be placed upstream of the first blade row a 
distance at least equal to two or three annulus 
heights. The best far-upstream inlet condition is 
straight axial flow with no wall curvature. Less 
iterations are usually needed for more widely spaced 
calculation stations; however, enough iterations 
should be used tc properly locate the streamlines. 
Calculation station spacing can vary somewhat along 
the annulus bui, as a general guideline, successive 
station increments should not be changed more than 
35 percent. 

When calculation stations are input close together, 
only some of them will be used for locating the 
streamlines if the station aspect ratio is above 7.0. 
This is done for program stability and convergence 
toward a solution. If the user does not specify which 
stations to eliminate from the streamline location 
procedure, the program has logic to do so when the 
station aspect ratio exceeds 7.0. The priority of 
stations kept for streamline location is as follows: (1) 
blade row exit stations are always used, (2) blade row 
inlet stations are kept if the blade row aspect ratio is 
less than 7.0, and (3) an annular station is kept if 
neither adjacent station is closer than the aspect ratio 
tolerance. 

The user can also specify that particular annular 
stations not be used for streamline definition through 
the alphanumeric station designation. The program 
looks for ROTO for rotor, STAT for stator, or 
ANNU for regular annular. Any other combinations 
of letters, numbers, or symbols designates the station 
.IS the extra-annular type. All the computations that 
are done for regular annular stations also are done 
for the extra-annular stations. The only difference is 
that the new streamline locations at that station are 
not used for the curve fit for streamline parameters. 
When the new curve fit streamlr.es are established, 
their intersections with the station line are found and 
the streamline parameters at that point are used in 
the equation -of-motion calculations. 

Ihe arrays of points that describe the hub and tip 
casing contours should extend at least from the 
furtlicst upstream calculation station to the furthest 
downstream one. There should be enough data 
points to adequately define the desired casing 
contours with a spline curve fit. 

The input boundary layer blockage factors have an 
option. A displacement thickness from the wall can 
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where 



be specified instead of blockage as a fraction of 
annulus height. This is done by using a negative 
number the magnitude of which is the value of 
displacement thickness. 

A total pressure profile can be input in place of 
losses. Although the way to activate this option has 
been discussed earlier, its full effects need to be 
understood. This option is activated for a particular 
blade row by using zero or a negative number in 
ILOSS (IROW). When the option is activated, an 
additional data card is required for that blade row 
(fig. 12(a)). The first parameter PTT(IROW), or P, 
in the equation, is the blade row tip f'arger radius) 
total pressure in psia. The five other parameters are 
polynomial constants P\ to P$\ therefore a total 
pressure at some other radial location is 

P = P,(1.0 + P,/? + P,R2 +P3/?3 


n -rn 

or the fraction of passage height at the blade row 
exit. Because these coefficients are stored into the 
locations of loss sets 4 and 5, those loss sets are 
destroyed for the run even if read in. 

When the pressure level is specified instead of 
losses for the last blade row of the compressor, there 
is an overspecification of data because the inlet 
pressure and compressor pressure ratio are input too. 
In computation the pressure ratio predominates; so 
tl.e pressure levels will be adjusted as necessary. Also 
note that when the pressure level is input, the total 
temperature profile must also be input (table II). 
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13 ► Siitiroutines use(^ r. input and iteration phases, 
figure Q. - Line representation o( subroutine calls. 


















Figured. - Concluded. 


At a rotor exit the total temperature level can be 
nnnt in place of the cumulative energy addition 

radial distribution of total temperature are input in 
th<> former oressure polynomial coefficient locations, 
pIKoW PRE(IROW), During regular 
Lation rhe program will um the polynomial form 


for rotor exit total temperature distribution when 
PRAnROWli> 100.0. The polynomial coefficient 
Stent^^^^^^ is found by adding or 

subtracting the number of IM's n^ded to give a 
remainder in the range - 100.0 to 1 00.0. 

When the total temperature level is input, 
pressure level can be set in two ways. It can either be 
determined from losses or innut directly by a 
polynomial, as discussed earlier in this section. 

The description of parameter variations i 
polynomials assures smoothness, but tne 
s^pecification of polynomial coefficients is not always 
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easy. In most cases the range of applicability for the 
polynomial independent variable is 0 to l.O. This 
considerably eases the burden on the user since 
computation is normally not needed to choose and 
set the polynomial coefficients. When the higher 
degree terms are used to define distributions, the end 
conditions are r ’atively easy to meet. However, 
some simple computations are needed to check the 
distribution. 

Another caution is that combinations of 
reasonable-looking numbers often give blade 
elements that one can judge to be poor by visual 
observation. The capability to make machine graphic 
plots of blade elements and the channel formed by 
adjacent blades is very useful. Such plots are made in 
subroutine EPLOT, which is activated by the input 
parameter OPM. Since graphics packages differ with 
computer systems, the program presented will not 
necessarily work directly on a user’s computer. 
However, i! is suggested that the user make the 
conversions necessary to plot the blade element 
surface arrays generated in EPLOT. 

The determination of acceptable polynomial 
coefficients for the centerline and thickness of an 
entire blade row can be difficult when high-degree 
terms are used. This task was eased considerably at 
NASA Lewis with an interactive graphics capability. 
A series of computer p-ogram.s were developed to 
design particular blade elements from actual 
centerline angle and thickness distributions. These 
data were then curve fit by least -squares methods to 
produce the input required by the program described 
in this report. Visual observaton of blade elements 


generated by this input for several fractions of 
annulus height is very helpful in avoiding obviously 
unacceptable configurations. 

The computer peripheral equipment also can be 
used by some other subroutines when options are 
activated with the pa^meter OPO. When the punch 
option is acti'aied, the tables cf fabrication 
coordinates shown on the listing are punched on 
cards in subroutine COORD, When the plot option 
of OPO is activated, subroutine BLUEPT plots 
tables of fabrication coordinates on a blueprint 
format. If a plot option is activated by either OPM or 
OPO, subroutine MERID is also called. It produces a 
meridional plane plot of the annulus flow path with 
the calculation stations and streamlines included. 

This code is interfaced with three other NASA 
codes through punched card output. Input for the 
TSONIC code (ref. 5), which is a blade-to-blade 
channel flow analysis code, is obtained with the T 
option of OPM. Input for the MERIDL code (ref. 6), 
which is a more detailed hub-to-shroud flow analysis 
code within a blade row, is obtained with the M 
option of OPO. Input for an off-design performance 
prediction code that is being developed at NASA 
Lewis is obtained with the O option of OPM. 

The computer program can be obtained from 
COSMIC, 112 Barrow Hall, University of Georgia, 
30601. The COSMIC program number is 
LEW-13505. 

Lewis Research Center 

National -Aeronautics and Space Administration 
Cleveland, Ohio, December 29, 1980 
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A 

annulus area; also streamtube channel area 

U 

local blade velocity, tt /sec 

'^1 

polynomial constants for as a function of 5 

u 

generalized variable in a dittereniiai 

a 

sonic velocity, ft/sec; also a coefficient in 


equation 


velocity gradient equation; also a 

V 

velocity, ft /sec 


polynomial coefficient 

V 

generalized variable in a dittereniiai 

b 

coefficient in velocity gradient equation; 


equation 


also a polynomial coefficient 

w 

weight flow, lb/ sec 

C 

constant 


a.xial distance, in. 

c, 

polynomial constants for conic radius as a 

a 

angle of streamline with reference to axial 

c 

function of S 

blade chord, in.; also a polynomial 

& 

direction, deg 

flow angle relative to meridional direction. 

Cpit) 

coefficient 

specific heal function for constant 

7 

deg 

blade chord angle, deg 


pressure, ft/sec^ 

6 

deviation angle, deg 

D 

blade element diffusion factor 

f 

angular coordinate on blade element layout 

D,.,= 

j simplified nomenclature, Dj= -(Q)/(/)/?^ 


cone, rad 

d 

r olynomial coefficient 

e 

circumferential direction, rad 

/ 

f.iciion force, ft /sec ^ 

K 

blade angle relative to local conic ray. deg 

H 

stagnation enthalpy, ft /sec ‘ 

X 

local angle of calculation station line with 

Up 

pressure-surface height, in. 


reference to radial direction, deg 
static density, slug/ft^ 
blade element .olidity, chord/tangential 

Hs 

suction-surface height, in. 

p 

0 

h 

static enthalpy, ft/sec- 


spacing 
time, sec 

i 

integer index; also incidence angle, deg 

T 

J 

integer index 

u} 

loss coefficient 

k 

curvature in curvilinear coordinate system, 
ft * ; also an integer index 

Subscripts: 

L 

distant'f along chord line, in. 

ca 

center of area 

1 

distance along calculation station line, in. 

/ 

calculation station index 

M 

Mach number 

i 

ideal value, as by an isentropic process 

m 

streamline direction in meridional plane. 

j 

streamline index 


in.; also an integer index 

le 

leading edge 

n 

streamline normal direction in meridional 

m 

streamline direction in meridional plane; 

P 

plane, in. 

stagnation pressure, Ib/ft" 

n 

also maximum thickness 
streamline normal direction in meridional 

P 

static pressure, Ib/ft^ 


plane 

R 

conic coordinate radius, in. 

o 

initial value 


series coefficients for polynomial, 

sp 

stacking point 


t 

transition point 

Rm 

radius of curvature in meridional plane, ft 

te 

trailing edge 

cR 

gas constant, ft lb/ slug 

e 

circumferential direction 

r 

radius from axis of rotation, in. 

1 

blade row inlet 

S 

blade element path distance, in. 

2 

blade row outlet 

s 

r 

entropy, ft sec- 
stagnation temperature, “R 

Superscript: 

t 

static temperature, *R; also blade elemen' 

( ) 

relative to rotor 


thickness, in. 

( )• 

flow at sonic condition (.V/' - 1.0) 


1 ? 


d 
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Appendix B 

Input Parameters for Compressor Design Program 

The input variables for the compressor design program and the associated options are 
appendix The format for the input data is given in figures 10 to 12. The calculation station and bla 
sets are input in the order in which they occur in the compressor flow. If any of the sets of option 
blade rows are needed they are considered part of the blade row set and they follow the particular b^ic blade 
riw dau set !n thTorder shown in figure 12 The only exception is any XCUT cards that are read in the output 
routines. These cards are at the end of the input data, but of course the sets of XCUT values must be placed in 

*^In folfowinVlist of param«ers^he independent variable S appears frequently. Since it is 
blade element definition variable, this preliminary explanation of its definition and . 

S in equations for the blade element centerline is the distance in either direction from the *^**‘®*^°" 
reference. The variable S in equations for the thickness distnbution is the dist^ce in cither - 

SSSm thickness point as a reference. All four of these usages of S are shown in figure 4. In all cases, S 
imshive aSTf^^^ their reference point. Tlie S values for thickness definition are normalized by 
blade elemrat chord. The S values for cemcrline definition are also normalized by blade 
mEF(lSw) ?s lis thii zero; however, when IDEF(IROW) is greater than zero. Sis normahzed to 1.0; that 

is, the maximum segment S is 1 .0. 
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(c) Stationary blade rows. 

Figure 11. - Input data format of calculation stations and basic blade row information. 


Parameter 


ACF(l.lROW). 

ACF(2.IROW). 

ACFO.IROW), 

ACF(4.IROW) 

ACRU.IROW). 

ACR(2,IROW). 

ACRO.IROW). 

ACR(4.lROW) 


Description Format 

This parameter is used twice to indicate options in alphanumeric form. As the A4 
first term of a data set it indicates the type of calculation station or blade row 
(ANNULAR, ROTOR, or STATOR). Any station description other than 
ANNU, ROTO, or ST AT will be treated as an extra-annular station, that is, th“ 
streamlines will not be forced to pass through the streamtube-fraction-of 
weight -flow point as determined by continuity at the station. The second use of 
AA later in the data set is the incidence angle option for blade design purposes. 
Inierpretable options are 2-D, 3-D, SUCTION, and TABLE. A 

noninterpretable incidence option word is set to the 2-D option. The 2-D and 
3-D options mean incidence angles are determined by procedures in reference 1 
for the respective option. The suction option gives zero incidence to the suction 
surface of the blade at the leading edge. The TABLE option means the blade 
incidence angles for the blade element will be input in tabular form. 
lNC(lROW,J). at the end of the data set. 

This parameter completes the incidence TABLE option discussed above. To A4 
reference incidence to the suction surface at the leading edge, the eight spaces ot 
the card for AA and AB must read 

TABLE SS 

AA AB 

(If ,AB is anything other than E SS, the incidence angles will be referenced to the 
leading-edge centerline.) 

polynomial coefficients for linear coefficient of blade element centerHne angle FlO.4 
equation for front segment, k = aS + bS‘ + cS^+ dS*^xh 
fl = ACFl -h ACF2 *R-hACF 3»/?‘ + ACF4*R-\ where R = (r, -r)/(r, -r^)— 
fraction of passage height at blade leading edge 

same as above for rear segment with same R FI 0.4 
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K L K < 3 , I R O W ) { F L F ( I , I H O W ) | F T K (1 , 1 R O W ) 


AT F ( 3. I R OW)|a T F (4, I R O W)|BTiF(l, 1 ROW) 


F T F ( 3, I R O W )Tc T F ( 4, JR OW)|p T F(l, I R O.W) 


AXK(4, I R QW)|B T R(I. 1 R 0|W) 


C T R (4, i R 0\V)|D T;R (1, I R O W) 


(cHfIDEF(IROW)>Q. 


, 1 R O \V)|d C F(3, 


B CR ( 2, I R O W)Ib C R (3, 


DCR4 2. 1 R OW)|d C R (3, 


. I R O W) E T E (3, 


B T F ( 2, I R O W)|b T F( 3. 


D T F( ;2, 1 R O \V)|dT F(3, 


B TR ( 2, • R O \V) p T R (3, 


I R O W)|D qp (4. 1 R O WA 


I R O W ) I B R (4 , I R 0 W ) 


I R O W ) I D q R (4 , 1 R O W > 


I R O W ) I E 'll E (4 . I R O W' ) 


1 R Q W)|b T|F (4, I R O W) 


I R O W) Id t|F(4, 1 ROW) 


! !r O W)|B 'nR(4, I R O W’) 


, I R O W)|D T|R {4. I R O W) 



1 M C (I HOW, 4 ) 


D E V ( I R O \V,3)|d E V (1 ROjW, 4 ) 


PHI { ;i R O W ,3)|P H I ( IR O'W. 4 ) 


TRANS OR O W,3) TRANSdROjW . 4 ) 


7. MAX (I ROW,3)|z MA X (I R,OW,4) 



if EE TABLE 


(d) If indicated parameters are TABU. 



X<* rjT|( 9 t I X FT T ( 1 0 » 


X F fH( 1 » 


(e) IfOPisVEL. DIA. 


F U T ( 4 4 1 X C F r ( 5 ) | X C I f ( 6 > 


X F V T i Nf''. 


X FI T - 4 1 1 X c i‘ r ( '■> 1 1 X 0 I T ( r> } 





(fllfNXCUTlIROWI< 0. 

Figure 12. - input data format of additional blade row information if needed by the options. 


Parameter 


Description 


Format 


ALIM(IROW) For a data set designated ROTOR, ALIM(IROW) is the minimum allowable 
relative flow angle (deg) leaving the rotor hub. For a data set designated 
STATOR, ALIM(IROW) is the maximum Mach number entering the stator at 
the hub. The program will reduce the stage energy addition to satisfy these 
conditions if a limit criterion has been reached during computation. If no 
aerodynamic limits have been reached in some other stages of a multistage 
compressor, the program will try to pick up the energy loss of the limiting stage 
in the stages free of aerodynamic limits. If all stages have reached some 
aerodynamic limit, the overall compressor pressure ratio is degraded to get all 
stages within the specified aerodynamic limits. The most efficient way to run the 
program is to specify the stage energy addition levels so than aerodynamic limits 
are not reached or at least not reached in a drastic fashion. 
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Parameter 


Description 


Format 



I 

r 


ATF(l.lROW), 

ATF(2.1ROW), 

ATF(3.1ROW). 

ATF(4.1ROW) 


ATRU.IROW), 

ATR(2,1R0W), 

ATRO.IROW), 

ATR(4,IROW) 

BB 


BCF(l,IROW), 

BCF(2,lROW). 

BCFO.IROW), 

BCF(4.1ROW) 

BCR(UROW), 

BCR{2,IROW). 

BCR(3,IROW), 

BCR(4,IROW) 

BH(1) 


polynomial coefficients for first coefficient a of blade element thickness equation F10.4 
forward of maximum thickness point 

£ = - «(' S. - S - ' - dS* 

with a = ATFl +ATF2»R + ATF3*/?- + ATF4*/?^, where R is fraction of 
passage height at blade leading edge and S„ is distance from maximum 
thickness point to centerline intersection of edge ellipse (fig. 4) 

same as above for reai i thickness with same R F10.4 


deviation angle option for blade design purposes. Interpretable options are 2-D, A4 
3-D, TABLE, CARTER, and MODIFY. Noninterpretable input is set to the 
2-D option. For the 2-D and 3-D options, deviation angles are determined by 
procedures of reference 1 for the corresponding option. The CARTER and 
MODIFY options are now the same in the program. They indicate the use of a 
Carter’s rule with a modification when the front and rear segments of a blade 
element have different camber rates I'he TABLE option means that the blade 
deviation angles for the blade elements will be input in tabular form, 
DEV(IROW.J), at the end of the data set. 

polynomial coefficients for quadratic coefficient of blade element centerline F10.4 
angle equation for front segment, k = k , + aS + bS' + cS- + dS* 
with b = BCFl + BCF2*/? + BCF3.R- + BCF4*/?’ , where 
R = (r,-r)/(r, -r/,)— fraction of passage height at blade leading edge 

same as above for rear segment with same R F10.4 


hub blockage factor for each calculaticm station; traction of the station annular F10.4 
area to be allowed for hub annular surface boundary layer blockage. The hub 
streamline will be displaced away from the physical wall a distance that gives the 
specified annular fraction. Negative input values are used as the magnitude of 
boundary layer displacement in inches. 


BM.ATLd ROTOR) rotor material density (Ib/in’). If a positive nonzero number is input, the blade F10.4 
will be stacked so as to balance out gas bending moniints with the ce.ntrifugal 
force moment for the material density. Because the hub stacking point stays 
fixed, the tip location is moved if necessary. 


BLADES(IROW) 

number 

BLEED(l) 

fraction 

BT(l) 

same as 


of blades in each rotor or stator blade row 

of weight flow bled off at particular calculation station 

BH(1) except applicable at tip 


F10.4 

F10.4 

F10.4 







w 







Parameter 


Descrip' ion 


Format 


BTF(UROW) 


BTR(l.lROW). 

BTR(2,IROW), 

BTR(3,lROW). 

BTR(4.IROW) 

CC 


CCF(UROW), 

CCF(2.1ROW). 

CCFO.IROW). 

CCF(4.1ROW) 

CCR(l.lROW). 

CCR(2.IROW), 

CCR(3,IROW). 

CCR(MROW) 


polynomial coefficients for quadratic coefficient of blade element thickness 
equation forward of maximum thickness point 


Zc 


ZC 


rai vSj, -S-vS„ + 


2v'S„ 


-bS' -cS^ -dS* 


with b = BTFl+BTF2*R + BTF3*/?- + BTF4*R3. where R is fraction of 
passage height at blade leading edge. 


F10.4 


same as above for rearward thickness with same R 


F10.4 


blade element geometry option for blade design purposes. Interpretable options A4 
are CIRCULAR, OPTIMUM, and TABLE. The CIRCULAR option gives 
circular arc blade elements. Noninterprttable input will be set to the 
CIRCULAR option. The OPTIMUM option means that the ratio of blade 
element segment turning rates will be set by an empirical function of inlet 
relative Mach number. Below an of 0.8 the blade element will be a circular 
arc. As is increased, the ratio of front segment turning rate to rear segment 
turning rate is reduced. A limit of zero camber on the suction surface of the 
front segment is approached at an Mj of about 1 .60. The TABLE option means 
the ratio of blade segment turning rates will be input in tabular form. 
PHl(lROW,J), at the end of the data set. 

polynomial coefficients for cubic coefficient of blade element centerline angle F10.4 
equation for front segment. k = k, +</S + bS^ +cS-’ +dS* 
with c = CCFl +CCF2»/? + CCF3*R- +CCF4-.«?\ where 
R = (r, -r)/(,r, fraction of passage height at blade leading edge 

same as above for rear segment with same R F10.4 


CHORDA(IROW). constants to define ratio of blade element chord to tip chord on projected F10.4 
CHORDB(IROW), plane 
CHORDC(IROW) ,. 

- - ^ 1 + /?*CHORDA(IROW) + R* *CHORDB(lROW) 

^ lip 

+ /?'*CHORDC{IROW) 

where R = {r,-r)/(r, -/•/,)— fraction of annulus height at blade stacking line 


CHOKE(IROW) desired minimum value of (A/A*)-\.0, where AM* is the ratio of local 
s/ireamtube area in the channel to the area required when M' = 1.0 within a 
blade passage. If zero is input, no adjustment will be attempted within the 
program. For input values greater than zero, incidence angle will be increased as 
necessary up to a maximum of + 2.0' on the leading edge of the suction surface 
in an attempt to give the specified choke margin at the covered channel entrance 
if the minimum occurs ai the channel inlet. 


CPCO(l) 
for 1 = 1 .6 


constants for specific heat polynomial function of temperature E20.8 

<7, = CPCOt 1 ) + CPCO(2)‘ T + CPCO(3)* l' + CPCO(4)* P 

+ CPC 0(5)» r* + CPCO(6)* P 


■1 









Parameter 


Description 


Format 


CRENGY 

(IROTOR) 


CTF(l,IROW), 

CTF(2,lROW), 

CTFO.IROW), 

CTF(4,IROW) 


CTR(l.IROW), 

CTR(2,IROW), 

CTR(3,IROW), 

CTR(4,IROW) 

DCF(l.lROW), 

DCF(2,IROW), 

DCF{3,IROW), 

DCF(4,IROW) 

DCR(I.IROW), 

DCR(2,lROW), 

DCR(3,lROW), 

DCR(4,IROW) 


DEV(IROW,J) 


DFTAB(K,J,1) 


DLIM(IROW) 


desired cumulative energy addition fraction through particular rotor to tot a. F10.4 
energy addition of compressor. Thus the fractions are progressively larger 
positive numbers through successive rotors. The last rotor must have 
CRENGY = 1.0 to meet the input pressure ratio. If a vaiue greater than z.u is 
input, the value is interpreted as a rotor e,<it total temperature level in degrees 
Rankine instead of the cumulative energy addition fraction. In the preexecution 
phase of computation the input temperature is converted and used as an 
appropriate energy addition value. 

polynomial coefficients for cubic coefficient of blade element thickness F10.4 

equation forward of maximum thickness point 

^ + -cs= 

with c = CTFl +CTF2*/? + CTF3«R2 ^.CTF4•/^^ where R is fraction of 
passage height at blade leading edge 

same as above for rearward thickness with same R 


polynomial coefficients for fourth degree coefficient of blade element cente.ane 
angle equation for front segment, k = k, +(^+bS +o +ao 
with = DCFl + DCF2*/? + DCF3*/?2 + dCF 4»/?3 . where 
Rs={r, —r)/{r, -r^)— fraction of passage height at blade leading edge 


same as above for rear segment with same R 


option control of location of transition point between segments of a blade 
element. The interpretable options l e CIRCULAR, SHOCK, and TABLE. 
The SHOCK option locates the transition point on the suction .suiiace at the 
normal shock impingement point from the leading edge of the adjacent blade. 
The TABLE option means the location of the transition point will be input in 
tabular form, TRANS (IROW,J), at the end of the data set. The CIRCULAR 
option and noninterpretable data put the transition point at midchord. 

deviation angle (deg) that can be specified by option. If the tabular option is 
used, a value is expected for each streamline stauing from the tip. 

blade element diffusion factor (D factor) for which profile losses are tabulated. 
Five values are input for each streamline; that is, K always has values from 1 to 
5, J is the streamline index, and 1 is the loss set index. The maximum number of 
sets is 5 Because D-factor values normally fall between 0.3 and 0.7, values of 
0.3, 0.4, 0.5, 0.6, and 0.7 for DFTAB on a streamline can be implied by leaving 
ihe DFTAB values blank. As a consequence of this option the DFTAB cannot 
be exactly 0.0 when K = 1 if you do not want the implied values of DFTAB. 

aerodynamic D-factor limit In a data set designated ROTOR this limit applies at 
the tip streamline. For a STATOR data set the limit applies at the hub. The 
program operates with this limit criterion in the same way as it did with 
ALIM(IROW) 
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Parameter 


Description 


Formal 




DLOS(KJ.l) 


DTF(l.lROW), 

DTF(2.!ROW), 

DTF(3,IROW), 

DTF(-l,lROW) 


DTR(l.lROW). 

DTR(2.1ROW), 

DTR(3,IROW). 

DTR(4.IROW) 

EB 


EE 


ELE(l.iROW), 

ELE(2.IROW), 

ELE(3,lROW). 

ELE(4.1ROW) 


ETE(I.IROW). 

ETE(2,IROW). 

ETE(3,IROW), 

ETE(4.lROW) 


Fl.OFRA(l) 


FFOVN(l) 


profile loss parameter lo cos /Ji 2a corresponding to DFTAB(K,J,1) reference 
arrays 

polynomial coefficient for fourth coefficient of blade element thickness equation 
forward of maximum thickness point 

L +J\s„-s-\ So + -4-r ) - bs- - cS-’ - ds-* 

V 2\5„/ 

with r/=DTFl +DTF2*/? + DTF3*/?’ +DTF4*/?-\ wheie R is fraction of 
passage height at blade leading edge 

same as above for rear segment with same R 


FB completes TABLE option of maximum thickness location. If the eight spaces 
controlling the option appear as 


T^EJJ 
EE EB 

the input values of ZMAX(IROWJ) will be used as the fraction of chord 
distance from the leading edge. If EB is not as shown, the values of 
ZMAX(IROWJ) will be used as the fraction of chord distance behind the 
transition point. 

option control of location of maximum thickness point of a blade element. The 
inierpretable options are TRAN and TABLE. The TRAN option and 
noninterpretable options will set the maximum thickness point at the transition 
point. TheT.^BLE option means the maximum thickness point location will be 
input in tabular form, ZMAX(IROWJ), at the end of the data set. 


coefficients for leading-edge ellipse . 

ratio of semimajor to semi minor ax^s 
minus 1 C 

=ELE1+ELE2*/? + ELE3*/?-+ELE4-/?-^ ^ 

a 

where R is ^Taction of passage height at blade leading edge 

coefficients for trailing-edge ellipse ratio of semimajor to semiminor axes minus I 

- -1 =ETEl+ETE2*/? + ETE3*/f2 + ETE4*/?' 
a 

where R is fraction of passage height at blade trailing edge 

cumulative weight-flow split between streamlines starting from tip. NTUBES, 
which is NSTRM -1, values are read. Thus the first value is greater than zero and 
succeeding values must increase to 1.0 in order for the last value to account for 
the accumulation of flow for all streamtubes. 

mass flow (lb sec) entering the firs; calculation station 



Jit 




F8.4 

F10.4 

F10.4 

A4 


A4 

F10.4 


F10.4 

F10.4 


Parameter 


Description 


Format 


IDEF(IROW) 


.LOSS(IROW) 


INC(IROW,J) 

MOLE 

NA 


NBROWS 

NHUB 

NLOSS 

NTIP 

NXCUT 


NSTRM 

OP 


blade definiticn index. When the index is zero, the blade segment centerline and 
surfaces are detined by d/c/dS = constant. When the index is not zero, the 
segment centerline and thickness are defined with fourth-degree functions of 
path distance from the transition and maximum thickness points, respectively. 

The specification of the coefficients for these functions is extra input, for which 
the format is shown in figure 12(c). If IDEF(IROW) is positive, the coefficients 
tor the definition polynomials are interpreted to be functions of segment length 
normalized to 1.0; but if IDEF(IROW) is negative, the coefficients are 
interpreted to be functions of segment length normalized by chord. The 
reference point for the centerline polynomials can be either the transition point 
or the segment ends. The possible combinations are shown in the IDEF(IROW) 
summary in table IV. 

designation of whicn profile loss set (I variable in DLOS(K,J,I)) to use with 15 
particular blade row. If the input value of ILOSS(IROW) is less than or equal to 
zero, a total pressure level is input in place of losses. The pressuie is input with 
the parameters shown in the first option of figure 12. These parameters are 
stored into the locations of loss sets 4 and 5; so those loss sets are not available 


for use with any blade row . 

incidence angle (deg) that can be input by option. If the tabular option is used, a F10.4 
value is expected foi each streamline starting frrm the tip. 

molecular weight of gas (28.97 for dry air) 15 

number of annular stations at which radial velocity profiles are constructed 15 
during computation 

number of blade rows (maximum of 20) 15 

number of points input to describe hub geometric boundary (maximum of 40) 
number of loss sets input (maximum of 5) 15 

number of points input to describe tip geometric boundary (m.aximum of 40). 
number of sections across blade for which fabrication coordinates are desired. If 1 10 


zero, the program will set the number of XCUT’s on the basis of aspect ratio. 

For all positive values the program will set appropriate locations to represent the 
blade. Negative values of NXCUT(IROW) trigger an option to read cards for 
the XCUT values. The number of values expected for a blade row is the absolute 
value of NXCUT(IROW). 

number of streamlines (maximum of 11) 15 

option controlling amount of output information desired. Interpreiable options A4 
are APPROX, VEl. DIA., DESIGN, and COORD. If the first four characters 
input in OP match none of the above, the program will try to proceed with the 
VEL. DIA. option. The program completes only velocity diagram information 
when run with the APPROX and VEL. DIA. options. With the APPROX 
option the locations of blade edges are estimated from the stacking line, but 
with the VEL. DIA. option the blade edge locations are input. The blade edge 
data are read from extra cards at the end of the data set for a particular blade 
type. The axiai coordinates are temporarily read into VTH(I,J), and the radial 
coordinates are temporarily read into PO(l,J). When run with the DESIGN and 


23 


Parameter 


Description 


Format 


OPM 


OPO 


PHKIROW.J) 


PRA(IROW), 

PRB(IROW), 

PRC(IROW), 

PRD(IROW), 

PRE(IROW) 


PO(I.J) 


PR 


COORD options, the program designs and stacks that particular blade row. 
With the DESIGN option only velocity diagram information is printed, but the 
blade leading- and trailing edge locations are for the stacked blade. The 
COORD option includes the printout of blade section properties and 
coordinates for fabrication. 


additional output options in effect if OP is DESIGN or COORD 


Card column 

Additional output 

[7 

8 

9 

10 


O 



Off-design punch 


T 



TSONIC punch 


M 



Blade element channel microfilm 


M 

O 


M and O options 


M 

T 


M and T options 


Additional output options in effect when OP is COORD 



Card column 

Additional output 

17 

18 


20 

1 

M 

i P 

! c 

M 

M 

1 

i 

j 

1 P 

! C 

i__ . - - 


Fabrication coordinate on microfilm 
Fabrication coordinate punch 
MERIDL punch I 

M and P options | 

M and C options J 


ratio of inlet segment turning to outlet segment turning (ratio of F10.4 

ydK/dS) 1 / (dK/dS) 2 ) for a blade element . If input values are expected by use of 
the tabular option, the data cards go with the optional cards at the end of the 
data set for each blade row. A value is expected for each streamline beginning 
from the tip. 

coefficients for polynomial equation to define profile behind blade row. Behind F10.4 
a rotor the pressure ratio profile is specified as 


— = 1 .0 PRA*/? r PRB»R‘ PRC*/?3 + PRD*/?-* + PRE*/?* 

■ ■ D 

where P, is the stagnation pressure at the rotor exit tip and a = 
(r, -r/,)— a fraction of passage height. When |PRA(IROW)l s IW.O, 
another option is activated. The input profile is for a temperature profile / // ^ 
instead of a pressure profile P/P/. The data value of PRA(IROW) is extracted 
from the input value by adding or subtracting lOO’s until the remainder is in the 
range of - 1(X).0 to 1(X).0. At a stationary biade row the polynomial is for the 
blade row exit tangential velocity profile in ft /sec. 

Pg = PRA//?2 -i-PRB//? + PRC-i-PRD*/?-(-PRE*/? 2 where /? = r/r, 


general stagnation pressure array in Ib'ft^ within program. The I 
station index and J is the streamline index. Only (PO(l,J), /"i. NSTRM) 
values are input; that is, the streamline value for the first calculation station. 
The input values are read in units of psia, 


When blade edge coordinates are input, some of the other PO(l,J) locations are 
used for temporary storage of the input values of radius. 


desired overall compressor pressure ratio 


F10.4 


F8.4 

F10.4 


I* 
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Parameter 


Description 


Format 


PTT(IROW), 

PTC(UROW), 

PTC(2,IROW), 

PTC(3,IROW), 

PTC(4,IROW), 

PTC(5,IROW) 

RHUB(I) 


ROT 

RTIP(I) 


SOLID(IROW) 

TALE(IROW), 

TBLE(IROW), 

TCLE(IROW), 

TDLE(IROW) 

TAMAX(IROW), 

TBMAX(IROW), 

TCMAX(IROW), 

TDMAX(IROW) 

TATE(IROW), 

TBTE(IROW), 

TC1E(IR07/), 

TDTE(:kOW) 


coefficients that describe blade row exit profile when it is input as an option. F10.4 
PIT is the blade row exit pressure in psia at the tip (highest radius). The other 
five values are polynomial coefficients for 

P= PTT»(1 .0 + PTCl *R + PTC2*/?2 + ptC3 + PTC4*R^ + PTC5-/?’ 
where R - (r, —r)/{r, —r/,) — fraction of passage height at blade row exit 

radius coordinates of a set of points that define geometric hub boundary FI 0.4 
(maximum of 40) 

compressor rotational speed, rpm FI 0.4 

radius coordinates of set of points that define geometric tip boundary (maximum F10.4 
of 40) 

tip solidity of a blade row (ratio of chord to circumferential spacing) FI 0.4 

polynomial coefficients of ra io of blade element leading-edge radius to chord, F10.4 

where rfe/c=TALE + TBLE*/J + TCLE*R2 + tDLE»R3 

where R (r^ -r*) — fraction of passage height at blade leading edge 


polynomial coefficients of ratio of blade element maximum thickness to chord, F10.4 
where r„a*/t^=TAMAX-HTBMAX*/? + TCMAX*R2 +tDMAX*R3 


polynomial coefficients of ratio of blade element trailing-edge radius to chord, F10.4 

where r,^/r=TATE-t-TBTE*/?-(-TCTE*/?2 +tCTE*R3 

where R {ji —r)/{ri — r/,) — fraction of passage height at blade trailing edge 


TILT(IROW) 


TITLE(I) 

TO(I,J) 


angle of stacking axis tilt (deg) in circumferential direction {r-6 plane). The angle 
is positive in the direction of rotor rotation. If |TILT(IROW)|> 100.0, a curved 
stacking line is specified according to r— /-„y^ = C(sin 7 -sin 7 ^^/), and the code 
of the TILT(IROW) is— 


0XXXX'XX 


tilt angle at tip in degrees. Circled digit 
controls sign of tip tilt angle. Even digit 
gives tip tilt angle same sign as hub tilt 
angle. Odd digit gives tip tilt angle 
opposite sign of hub tilt angle. 

For example: 12332.65 gives a hub angle of 23° and a tip angle c' - 32.65°. 

description of compressor for printout and later identification 18A4 

general stagnation temperature array in program. Only (TO(l,J), J = l, F10.4 
NSTRM) values are input; that is the streamline value for the first calculation 
station. The input values are in units of °R. 



tilt angle at hub in degrees with sign of 
overall TILT(IROW) number 
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l\ir;imctov 


Dcscnpp.on 


K>!i1UU 



rRANS(lROW..n locatu>n of tiansition pomi i>n blade demcni centerline as traction ot blade I 10.4 
element chord. It input wiliies arc expected by use ot the tabular option, the 
data cards gi> with the optional cards at i!ie end ot the data set tor each blade 
row. .A value is expected tor each streamline beginninii trom the tip. 

VIH(I.J) general tangentiai component ot velocity array in program. Only (VlH^l.J). MO. 4 

J 1. NSTRN!) values are input; that is, the streamline value tor the tirst 
calculation station. The input values have units v>t‘ tt sec. 


When blade edge coordinates are inpui. some of the other VTMvI.J) locations F8.4 
are used for temporary storage of the axial coordinates of the pciints. 

XC'UT(K’) radial location of blade section planes. Whether or not data cards are read for FM0.4 

values of XC^U T(1C) for a blade row is controiled by the value of NXCllT {lO. 

Any XC11T(!C'; cards are read in an output routine. Therefore they must follow 
all cards read in .subroutine INPUT; that is. they follow the ANNULAR card 
for the la.st calculation station. There is no index identifying the data with a 
particular blade row, so the data sets for the blade rows are expected in the 
order that one would sec the blade rows in moving through the compressor from 
the inlet. Start the set of points for each blade row on a new card. It is 
preferable, but not necessary, to list the XC'UKIO lor a blade row in order 
starting from the tip. 

XHl^BlD axial coordinates of set of points ihat define getvmciric hub boundary. The axial FIO.4 

extent of the coin dinates must ai least reach the first end last calculation 
stations. The hub coordinates must have the same rrfcrcucc origin as other 
input axial coordinates, thai vs. casing, olade edge, and stacking line 
coordinates. The number of points input should be 4<.fi s4?> 

XTlP(l) axial coordinates of set of points that define geometric tip boundary (See I 10.4 

XHUB(I) fv>r additional comments.) 


/MUH(l) blade data set hub-axial coordinate. When the data set is a blade rather than an MO. 4 

ANNUl AR station. /HU B(l) is the axial location of the blade stacking line at 
the hub. 


/M AX(lROW .J) location of maximum thickness point as fraction of blade element chord. If input F10.4 
values are expected by use of the tabular options, the data cards go with the 
optional cards at the end of the data set for each blade row. A value is expected 
for each streamline beginning from the tip with a leading*edge or transition- 
point reference according to option (see FB). With a transition point reference 
the values input are (m - /) c 



/TIP 


blade data set tip axial coordinate. (Sec /Hl’B(l) for similar additional F10.4 
comments.) 
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Appendix C 

Development of Equations of Motion into Form 
Used in Computer Program 


In itic coinpulcr prt>(irain ihc equations ot motion 
are applied at ealeulation stations that arc presumed 
to be outside the blade rows; so tlie equatiotts ol 
motion are more conveniently developed in an 
absolute, rather than a relative, coordinate system. 
The general equation of motion (eq. .M21) of ref. 7) is 


In tertns of equatiotis and (C'6) the term 

r V, ( V \ I ) can be expressed as 


I xt V X n = 


() rh ri 


+ V// = f x(V X n + rv.v+/ (01) 

()r 




m 

dn 


■b 0-A', 


rdn 


When steady llow is assumed and the local friction 
force is ignored, equation (C'l) reduces to 


V//= fx(V X f*) f fVs (O:) 

In orthogonal curvilinear coordinates the velocity 
vector can be expressed as 


.. (»10| f rh|Ol + 


r dn 




dn 


^ f 


(C7) 


f - r-dlf (<- .)) 

where m is in the streamline direction in the 
meridional plane and n is in the normal direction in 
the meridional plane. Of course f „ is zero 
ywhere for this application. I he curl term in 
general can be expressed as 


Now break equation (C2) into the three component 
equations. In the ^ direction 


d// 



(08) 


The zero in equation (C 8) recognizes circuintcrential 
symmetry of .v. In the meridional plane streamline 
direct ion 


d(rf g) 

df^, 

n 

A, 


_ d(r) ^) 

3n 

dd 

T 

r 

.W 

dm 


wkeie k„ and k„ are the curvature of the streamline 
and the normal, respectively. All terms containing 
y„ are zero for this application. The a.ssumption ol 
symmetric flow in the circumferential direction 
makes .If „,/(W equal to zero, .'kiso, because angular 
momentum does not change on streamlittes outside 
the blade rows 


<1// _ Av 
dm dm 


(t*)) 


The zero in equation (0*)) comes from the 
assumption that entropy does not change along 
streamlines that are outside the blade rows. In the 
meridional plane normal direction 


dH 

dn 




dn dn 




dn 


(CIO) 


(t) 

3m 

Thus equation (C4) reduces to 


V X 



dn 




fh d(rf(() 
r dn 


Fquations (C8) to (CIO) apply to the three curvilinear 
component directions. However, in the program 
velocity and state values arc available along station 
lines; so it is of computational convenience to apply a 
component equation along a station line. To 
iiccomplish this objective, the derivatives in Ihc 
meridional plane are converted from the orthogonal 


p- 



r 


streamline aiul nonual dncctions to the gcncrallv 
noiUMihOfcioiial strcainluie ami staium liiie Jireetions. 
I he an>ile nomenelature tor the conversion is shown 
m tijiiire l .V 

1 he cnihalpv iiiavlieiu m the statu n line dirc\tion 
can he evpressevl as 


1 heretoie 

til cos(tv ^ \) 


^ HI 

" ^ cos(t^ f \) ♦ * sm(i^ f \) 

,7/ iVr 


.Hi iit^ iHi iim .Hi Jn 
r iki iil ^ ,krt ,il ^ ,// 


:0| ;0i ♦ |0i'.sm(iv ♦ \) ♦ * cos(«» <■ X) 
tin 


</// .W 
<// " .1« ' 


eos(tv f X) 


rhcrct'orc 


in’ dy 


di cosiiv ^ X) dm 


tan(n > X) {C\4) 


ds dv . ^ . vv 

. cos(o ♦ X) f _ sin(o! > X) 
dl dn dm 


In general a stativ>n line derivative can he expressed as 

d _ dn d dm 
di dn di ^ dm di 

vv 

^ cos(tt f X) f sm(u f X) [CM) 

dn dm 


^ . cos(n + X) f |0|sin(<« + X) 

dn 


Fheret'ore 


dx _ ds \ 
dn di cv>s(i^ r X) 


When equation (C'12) is applied to the other normal 
dc»i\aiivcs of equation (CMO), the lollowing relation 
de\ clops: 

c/(rl _ d{rl ft) d(rl fi) 

,, V cos(tv t X> ♦ sin(ii f X) 

di dn dm 


The application ol equations (Cl 2) thn>ugh (Cl 5) 
to (C IO) gives 


j! ^ u ^ ^ .Ji ^ ff* stn(«n f X) 


di r di 


d{r\ ft) 


;os(tt f X) ♦ 01sin(cr ♦ X) 


I • A ,„cos(it f X) t r 


The stieamhnc curvature is 



Mjltion I Mir’ 


Axi.ll JlMStliVl 




, iXk I -,-^v 

where R,ft meridional plane streamline radius 
of curvature. Substituting equation (C'I7) into (C'I6) 
yields the following form for the meridional veUvity 
gradient: 

I m _ dH d(r\ dl ^ 

dl - dl rdl dm 

♦ o'” cox(.» ♦ \) (C'18) 

r\ (// 


U\c state pFx j»vTties api>earing in equation (C'I8) are 
//, r, and a\ However, two st^te properties are 
sufficient to establish the others at a pivint. For a 


thermally perfect gas (p=(9tRt) >t is rather e^y to 
compute other state properties from two selected 
properties; . o it is desirable from a computer storage 
standpoint to store only two properties throughout 
ihe flow field. The two properties selected were 
stagnation temperature and pressure. These two 
properties, along with the velocity components, are 
sufficient information for the calculation of the 
state properties. If these two properties can 
directly in the equations of motion, the need to 
compute some state properties may not exist. To 
express s in terms of f and P. start with the property 
relations 

P 

For the introduction of stagnation properties note 
that the thermodynamic process of moving between 
the static and stagnation states is isentropic by 
definition. Thus equation (C19) for this process 
becomes 


(is dh _ I 
til ^ t dl ~ pt 


dj^Pe.t/r| - l/^w(^^ .>(0//dt 


^ _ I exol - i 

di t dl pt dl [ tR 'r l 




1 dh i rfP/ PU (^^^dtl 

" i dl pt di^P’ (Sifit\P'dl[h t J 


1 ^ _^dP d 
1 dl P dl dl 


Jt t J 


For a calorically nonperfect gas this becomes 


The application of Liebnitz's rule to the last term 


= cAt)d( 


p f 

J I) p J f t 


> 1 

'■ iR )/ t 




Equation (C19) used as a derivative with path 
distance can be written as 




C (T) dT _ ^ dt 
T dl t dl 


The variable fc„U)/r) is not a direct function of path 
distance; it is a function of temperature alone. 
Therefore the partial derivative with respect to 
distance must be zero. Thus the derivative of the 
integral can be expressed in terms of gradients at the 
limits so that 

di[]( t^’r T dl t dl 


Substituting equation (C20) gives 


Substituting (C23) into (C22) gives 








iis 

lil 


I dh 
t dl 


dP 1 dH 
P d^ f dl 


1 dll 
i dl 


ds ^ I dH 
dl ~ T di 


dP 
P dl 


I dH 
f dl 


1 dP 
di 


(C24) 


Equation (C24) is essentially equation (C21) 
expressed in stagnation state variables. Equation 
(C24) would turn out to be the .same for a calorically 
perfect gas. Substituting equation (C24) into (CIS) 
gives 


dV„ 


dl 


1 / dirk'^) . 


dl 


yi 


+ cos(a + X) - - --- + 


T dl P dl 


A rearrangement with all the state property terms 
together gives 




dl 




dl rdl 


dl 


9V„ 




(C25) 


•All the terms on the right side of equation (C25) 
can be computed quite accurately except dV„/dm, 
which is the gradient of along a streamline in the 
meridional plane. The distance over which dV„/dm 
changes sign are of the order of the calculation 
station spacing so that representative values of 
dV„/dm cannot be obtained from a V,,, curve fit 
along meridional .streamlines. A better value of this 
derivative probably can be obtained by means of 
local continuity. From equation 9(12) of reference 7 
differential continuity can be expressed as 


I dp 


p Dl 


+ ^r=o 


However. 


I Dp 
P Dl 


1 Dh 
a' Di 


.10 


so equation (C26) can be written as 
1 Dh 


Di 


+ VK=0 


(C27) 


Equation (C27) expanded from its vector form is 


\dt r 


'"am ^ '•any 


^ 1 d{r V„) \9y» \ a^) 

r dm r dd ^ r dn 


4 " ^ m^m 4 " ^ n^n 


Outside the blade rows the flow is assumed to be 
axisymmetric and steady. Also, because there is no 
velocity component normal to the streamline, the 
equation reduces to 


dm r dm 


(C28) 


Stagnation enthalpy is defined as 


H = h+ 


(C29) 


dH dh 


_ »,• , y ^0 

dm dm dm dm 


But because a///am = 0 outside the blade rows. 


= _ I.' _ f 

dm "'dm * dm 


(C30) 


Outside the blade rows angular momentum is 
conserved along streamlines; sc 


(C26) 0 = + r 


dm dm ” ' dm 
Rearrangement gives 


dVe 


y» dr 




— — ~ — — sin « 

dm r dm r 


(C3I) 


Substituting equation (C31) into (C30) gives 


Substituting equation (C34) into (C 33) gives 


■ = -K„ 


dm r 


lLul = 

dm 


M - - IL '* 


Substituting equation (C32) into (CIS) gives 

— . - Cm + — Sin a ) + — ^ 
g2 \ ’” dm r J r dm 


+ ^+^'«*m=0 


(l - sin rt+ V„k„=0 

\ } dftX V ^ 


9m m\, 


_ sin a + K^*„] (C33) 


The curvature of the streamline normal k^y which is 
dot/driy needs to be expressed in terms that can be 
evaluated. 


^ ^ cos(a + X) + ^ sin(a + X) 

dl dn dm 


da da I da sin (a ^ X) 

^ ~ ^ cos(aVX) dm .'Os(a + X) 


+ ^ sec(a + X) 
dl 


tan(a + X) 


Calculation of dVf^fdm by using equation (C35) 
should give a somewhat more accurate result than a 
curve fit or a finite difference computation across 
increments that span whole blade elements. 
However, a potential divide-by-zero complication 
has been introduced with the term — 1 . In 
equation (C35) the term in braces in essence 
represents the dA/A term of one-dimensional now 
theory. At a Mach number of l.O, dA/A is zero, 
which is the throat of a nozzle. For compressor blade 
rows the throat occurs within the blade passages. 
Internal flows adjust around locally choked regions 
so that the throughflow Mach number outside the 
blade only approaches 1 . Computation of the 
detailed nature of the flow is not available from only 
stations outside the blade row; so a minimum value is 
imposed on the denominator through an empirical 
additive term to help siabilize the iterative procedure. 
The additive center term is 




/=0.1^-expj-10(M^-n 




da da , , ,, tan(a + \) 


Its characteristics and effect on the denominator are 
shown in table V. 
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Appendix D 

Conic Coordinates of Blade Centerline Path 


Local blade angle is defined with respect to the 
local conic ray (fig. 14). Let the blade angle vary with 
path distance along the cone according to the 
polynomial 

K~Kf 4-c5-^ -hdS** (DI) 

where kj is the blade angle at the transition point 
between segments in this application. The path 
distance S is with respect to the transition point 
reference but always positive in the direction from 
inlet to outlet. 

The conic radial component of the centerline can 
be found by integrating the differential equation for 
that component 



Figure 14. - Blade element centerline nomenclature. 


CCS #C = 1 ~ 


2! 




(D3) 


rf/? = cos|#cldS=cos(K, -haS+frS^ -PcS^ ^dS^)dS 

(D2) 

The problem is that a trigonometric function of a 
polynomial is not readily imegratable in closed form. 
However, the ^unction can be expanded in series 
form and integrated term by term. Of course the 
series is infinite but it is convergent within the range 
of our application. In the following presentation 
enough development is given to show the form of the 
series. Upon application in the program a tolerance is 
used so that no more terms than necessary are 
calculated. 


When equation (Dl) is substituted, the terms of like 
powers of S can be summed to give in symbolic form 

cosx = | || +1 IzS^l I3S-+I \4S^ + . . . (D4) 

cos K ds = I !|S4"j b ^ 

+ i +1 14 j--- 

When terms of similar coefficients are combined, the 
following form evolves: 


s 


cos K dS = - 
a 


cos Kj sin(aS) 4- sin Kf cos(aS) 


1 . 

— sin Kf 
a 


I 


1 


n 



-hbsin 
+ 6 cos 
h- 

+ y COS Kf 

ft* . 

+ sm K, 


^ ^ S5 _ S’ S’ 

T T y “ ^ T 6!' T 

S* S* j* S* \ 
“T IT y ~ 5! T‘ • V 


\ 5 2 7' 4! 9 ■ ■ 7 



^5! ^ 

3! 8 ■ ■ 7 


■) 






j 

i 

I 

i 


I 

it 

i 


/>•» 

+ — r COS K 

4! 


/)•' . / S' a~ \ 

+ ^ - y -^...) 

ft’ / S« \ 

ft-» (S'^ \ 

+ -jC0SK,(^ - ...; 

+ ftcosA,[-c-+^ ~+(^ --— j — 

a^c a-^c^ [ a^c a^c^ls*^ 

L ~ JuT) ''■ IT J IT L “ 'sT ■'■ ^ “ 2 ( 3 !) J 1 ? 

^ . ( S^ /c2 ah\s^ r flV a’cls" 

ft- . / S« a^cS"’ \ 

"T*'"'‘'V"T-T- lo ^■• 7 


ft- [ sw a’c- c^\ 

yCOS.,[.C^+^-~ + -j 

ft’ / S'« \ 

+ -jcos.,(^.-ro + •••/' 


^\S" 


ft’ s" 

JT sinA,|^-tfc- y- + 


^ [ _,S' fl-JS*’ ad- S'= a^dS'-’ 

+ ftCOS«,[-dy + -y ^ - y- -- + - -y + — -- 


o’d2 s'-* / a»d d’\s” 


3!{2) 14 8 ■'■ 3! / 15 


. . ,5" a’dS"’ d^ S" a’d S'2 a^d- 

+ ftsin«, ad g 3, 10 ^ 2 11 5! 12 2(2) 13 


_ a ' s'-* a_-*d- S'^ l" ^ V 

7! I 4 41(2) 15 '*’\ 9! 3! / 
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. /'S’ fl2^S" a*d fl</2 J14 

+ -- ^ - 


2 

b^ 


) 


‘•'G 


9 2 11 '' 4! 13 2 14 6! 15 ^ 

S»0 flVS>‘ , rf2 513 ^ ^ \ 

! 14 2 (2) 15 ■ 7 


+ — cosK,lad 12 + 2 13 5! 14 2 (2) 


/ S" 

-COS.,(^rf-jy- 


a^d S'^ a*d S*^ ad^ S'* 


+ 


2 13 4! 15 


Si6 \ 

16 ■ 7 


b^ . ( jS‘2 a^d Si4 S15 \ 

3T"‘n«rl^-«rf-i2 + TT 14 " T TT ^ • 7 




- cos*,\^-ad-jj-t-. ..j 
b^ ( .S'5 '\ 

5T'°”'V"'T5"" ■■ ■) 


+ c sin #c 


( S* a 
'\ 4 ■'■ 2 


2 6 4!' '8 


/ S5 a3 fl* S’ 

+ c"COSK^^ ^ ^ ^ y *7 51 

a?- \ 

+ yCOSK,l^-y + y ...; 

+ ysin«,(^ay+...j 

( 4 . 


+ C cos X, 


S« . S^ g</^ S'3 

10 ''' 4! 12 2 13 


•) 


+ c Sin X 


/ ,S’ flJdS“ rf2 S>2 ^ flV S>3 \ 

'V^ 9 3! 11 ^ 2 12 5! 13 ■ 7 


+ 2 sinx,|^d jj 2 13 ■ 7 
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y cos«,(ad^+. . .) 


+ </sinKj^ j 2 7 4! 9 


. S6 c3 S8 ^ 

+ dC0SK,l^-fl-g+ 3| y + 


3! 8 / 


rf2 

+ — COS K/ 


(-?•) 


Cs“ a2 S'3 ab S^* ^ 

+abcdcosK,^jY 3; 13 2(2) 14 

r a^b _ ^1 ^ 
■^L4!(2) 4 J 16 


0^ 

IT "3K4)Jir 




a^c 


7! 4!(2) (3!)2 


g2f>2 c2]5»2] 

/3n2 3! J 17 j 


f aS'2 &S'3 c\S‘4 

■abcdsinx,l- 2/14 


o2ft 

3!(2) 15 
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fl&2 


6! 2(3!) 3 


r a^b abc b^ 1 ] 
)!(2)JT6"^L~5!C)''’ 8 4'Jn) 


d- ( 

' S'5 fl2 


+ abc y sin 

. 15 3! 

17 


/ a S>6 


-^abc^ cos Kf 

\ ■ 2 16 

...j 




With these groupings shown, patterns of terms and 
coefficients can be observed. The whole equation was 
coded into three rather brief subroutines — one for 
terms with two coefficients, COEFl (two of the four 
coefficients a, b, c, and d); another for terms with 
three coefficients, COEF2; and one for terms with all 
;.,"ir coefficients, COEF3. Finally the coefficients of 
the terms with the same powers of S are summed: so 
the [ 1 terms are known in 
S2 

R=R, + [ ),S + ( )2 y +( h J 


4(14^+.. -MUf 

Because in the following developments these 
coefficients appear frequently within parentheses, 
for simplicity the [ ]’s are replaced with c’s; that is, 

52 53 54 S" 

/? = /?, +rjS + C2 y +C3 y +C4 — 

yD6) 
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The conic angular cocdinate can be expressed as 





^ sin K 
0 ~R~ 


dS 


(D7) 


where both sin k and R can be expressed as infinite, 
but convergent for our purposes, polynomials of S. 
Since a polynomial in the denominator is an 
undesirable form to integrate, the polynomial for R 
was converted to a polynomieil in the numerator of 
the form shown in equation (D8). 


The conversion from equation (D6) to (D8) begins 
as 


^ ^ 

^ ^a+CiS+C2(S2/2) + C5(SV3) + . . . 


J 

1 +(C|//?,)S + (C2//?,)S2+(c3//?,)S3 + . . . 




sin K 

0 “iF 


dS 


= 1 

R. Jo R 


sin K dS 


where 


^ 1 

1-Z)/S-Z)252-Z)3S3_ 

where 


R 


= 1+/?,S + /?2S2+/?3S^+. . . 


(D8) 




I -DiS-DiS^-D^S^ 


Rl /?3 R^ 

1 ^DjS + iDi +Of)S^ ■*- (£>3 2D| Z)j +(£>4 +2£)|Dj +£>§ + 3£)|D2 + 0|)S^ 

! 

^ I -PiS-DjS^ -Z?j53 

OiS + D2S^ +D}S^ +£>4S^ 

D,S-DjS^ -D,D2S^ -D,0,S* 

(D 2 + Df)S2 + 0|Z?2)S3 +(£>4 +£)|£)3)S< 

(£)j + 2£)|Z)2 +£>|)S3 +(£>4 +D|£>3 +£)|+D2Df)S^ 

(D] +2£»i£> 2 +Z)f)S3 -£>|(£)3 +2£>3 +2£>i02 +Df)S< 

(Di + 2D, D3 + D^ + 3DJD2 + D^ )S^ 
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Table VI summarizes the preceding division. 

The coefficients for equation (D8) are generated in 
subroutine RCOEF. The coding for the procedure is 
somewhat complex, but in general not much 
computation is required to satisfy a tolerance 
criterion of l.OE-08. 

The conversion of sin k, where 

k = k,+cS + 6S2+cS3 +^5^* 

to the polynomial form 


1 fS/?, . 

= lj^(l+/?/S + /?2S2+/?jS3 + ...) 

X iA -f- >1 2 “f” 3 4 + . . . ) 

^ l)S 


sin x = -f > 445 ^ . . (E)9) 

is accomplished in the same way as it was for the 
cosine series (eqs. (Dl) to (D5)). In fact, the cosine 
series can be converted to the sine scries with the 
following substitutions: 


Cosine scries 

-sin Kf 
- cos K( 
sin Kf 
cos K; 


Sine scries 

cos K( 

- sin K( 

- cos K( 
sin Kf 


■f (.<44 '¥R\A'^ R2A2 R^A I )S^ 


, -43 +/J;>42 +/?2-4i J.J 

4 ~ i 


4 ^4 +^1-^3 +^2^2 +^3^1 ^4^. 




Consequently the same routines that are used to 
compute the cosine series ca^i e'lsily be modified to 
compute the sine series coefficients also. 

When the polynomial series coefficients in 
equations (D8) and (D9) are known, the integration 
for e is straightforward. 


The general routine for establishing the 
polynomial coefficients for the conic coordinates is 
EPSL2. The end result is constant polynomial 
coefficients for the conic coordinates (/? and €) as a 
function of S. These coefficients are saved so that the 
conic coordinate at any S of interest can be computed 
easily with subroutine CONE. 
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T.4BLE I. - OVERVIEN^' OF COMPUTER PROGRAM 



Program control 


Input ar.d initializaiion 

Iteration 

Terminal calculations 

Head and interpret data 

Outer loop: 

Overall blade row performance 

Locate calculation stations 

At calculation stations 

on streamlines at calculation 

Set coefficients of equation 


station: 

At each station for ea?h streamline, 

of moUc' 

General 

estimate stagnation temperature 

If blade design option, set 

State properties (temper- 

and pressure and axial and tan- 

incidence nnd de\ iauon 

ature and pressure) 

gential velocities 

angles, compute new blade 

Velocity diagrams 


edge locution, and reset 

Streamline information 


calculation station location 

1 

1 

Inner loop: 

At each calculation station 
Solve for meridional veloc- 

Blade rows 

Element definition param- 
eters 

Incidence and deviation 
angles 

Aerodynamic performance 
parameters 


ity distribution to satisfy 
equations of motior and 


continuity 

Reset streamline locaUon 

Streamline choke margin 

Blade section parameters; 
Surface coordinates 
Area, moments, etc. 


I 
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utnUal velocity conriponent at stator exit 



St<^na»ion pressure at tip (psia) - PPJflHtAV) 



TABLK in. - KXAMPLK PKOBLKM 




L 


4 ! 


PROFILE LOSS TABLE NO. 





TIP AXIAL LOCATION HUB AXIAL LOCATION TIP BLOCKAGE FACTOR HUB BLOCKAGE FACTOR MASS BLEED FRACTION 

(IMCHE'> riMCHES) 


TABLE III. - Continued, 




4.1 


I 


■ 



TABLE III. - Continued. 
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PPINTOLIT OF IMPUT STATION DATA *** 




TABLfc III. - Continued. 



0000 9 3000 1.0000 0.0000 0.5000 
0000 9.6u00 I. 0000 0.0000 0.5000 
QOOO 11.1000 1.0000 0.0000 0.5000 
0000 16.0000 1.0000 0.0000 0.5000 


TABLK III, •“ Continued. 
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